Universality of attractors at weak dissipation and particles distribution in turbulence 
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We study stationary solutions to the continuity equation for weakly compressible flows. These 
describe non-equilibrium steady states of weakly dissipative dynamical systems. Compressibility 
is a singular perturbation that changes the steady state density from a constant " microcanonical" 
distribution into a singular multifractal measure supported on the "strange attractor" . We introduce 
a representation of the latter and show that the space-averaged properties are described universally 
by a log-normal distribution determined by a single structure function. The spectrum of fractal 
dimensions is derived. Application to the problem of distribution of particles in turbulence gives 
testable predictions for real turbulence and stresses the role of pressure fluctuations. 
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Phase space density, that describes averages of dynam- 
ical variables, is one of the most fundamental objects of 
the classical equilibrium statistical mechanics 0, [|[ . It 
obeys the continuity or Liouville equation which time- 
independent solution describes a steady state. For mix- 
ing systems the solution is the microcanonical uniform 
distribution. A constant solution exists due to the in- 
comprcssibility of the velocity field in the phase space. 
Here we study a perturbation of the above scheme in the 
general setting of weakly dissipative dynamical systems 
[2, Q • A dynamical system is represented by a point in 
a d— dimensional space which motion is governed by the 
local value of a prescribed smooth velocity field. Weak 
dissipation corresponds to small compressibility of the 
velocity. We allow for a generic, not necessarily Hamil- 
tonian, incompressible component and the results apply 
both to a time-independent deterministic flow and to a 
time-dependent flow which statistics defined by spatial 
averaging is stationary. 

Compressible perturbation is singular. It transforms 
the stationary density from a constant to a singular den- 
sity, traditionally called the Sinai-Ruelle-Bowen (SRB) 
measure supported on a multi-fractal. This needs a 

special approach which is developed here. As in the usual 
near-equilibrium physics, the properties of the SRB mea- 
sure at weak compressibility are expressed via the fluc- 
tuations in the unperturbed microcanonical ensemble. 

The states described by the SRB measure have a finite 
rate of the Gibbs entropy production. The rate finite- 
ness does not contradict that in the steady state the 
entropy should be constant, as the latter is not finite 
by the measure singularity. Thus the SRB states are a 
model of non-equilibrium steady states where the system 
exchanges entropy with the environment constantly. 

An example of a dissipative dynamical system that 
played a crucial role in the development of understand- 
ing of chaos and non-equilibrium physics is the Lorenz 
oscillator [6]. It demonstrated what is now known as a 
general property of the dissipative dynamics. At large 
times the trajectories asymptote a zero-volume, multi- 
fractal set in space - the "strange attractor". Thus a 
small perturbation transforming incompressible velocity 
into a compressible one leads to a non-small change in 



the steady state density: it becomes infinite on its multi- 
fractal support. Such a dramatic change is due to the 
infinite time of evolution implied by the consideration of 
a steady state. While locally compressibility leads to a 
small disbalance of trajectories going in and out of phase 
space regions, the disbalance effect accumulated over a 
long time is big. The analysis should therefore disentan- 
gle the effects of long time of evolution and of small local 
disbalances that can be tackled as a small perturbation. 

For a mixing incompressible velocity the evolution of a 
small volume of the phase space makes it dense in space. 
Its coarse-graining over an arbitrarily small scale covers 
all the available phase space, assumed finite here. When 
a small compressible component is added to the velocity, 
the coarse-graining of the evolved volume over an arbi- 
trarily small scale does not cover the whole space any 
longer. However the coarse-graining over a small but fi- 
nite scale, that tends to zero analytically with compress- 
ibility, already covers the whole volume. Thus at small 
compressibility the attractor is "almost" dense in space. 
The gaps in its structure have volume equal to the whole 
volume of the phase space but a minute coarse-graining 
is enough to fill the space. It should be stressed the at- 
tractor does have a non-trivial structure. 

A practically important example of weakly dissipative 
dynamics is heavy particles in a turbulent flow. Here 
physical and phase spaces coincide and one can actually 
see the attractor. The latter, like the velocity, evolves in 
time keeping its space- averaged properties constant. At 
small compressibility, universality allows to describe the 
particles distribution in real turbulence without knowing 
the details of the statistics of turbulence itself. 

We pass to the analysis of the continuity equation 



d t n + V • (nv) = 0, 



(1) 



for the phase-space density n of a general dynamical sys- 
tem which trajectories q(t, x) are governed by a smooth 
d— dimensional compressible velocity field v in a finite 
domain (which volume we set to unity) according to 



<9 t q(i,x) = v[t, q(t,x)], q(0,x)=x. 



(2) 



We study both the case of a time-independent determin- 
istic flow v(t, x) = v(x) and the case of a time-dependent 
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v(t, x) stationary with respect to the statistics defined by 
the spatial averaging. The steady state solution Usrb 
to Eq. ([IJ is time-independent in the former and time- 
dependent and stationary in the latter cases. We study 
the space averaged characteristics of tisrb , that are time- 
independent in both cases. The analysis is written for 
v = v(x) while generalizations to v = v(t,x) case are 
either straightforward or follow 7, 8] literally. We as- 
sume v is weakly compressible and its decomposition into 
transversal u and longitudinal u' components obeys 

v = u+u', V • v = V • u'=u;, v! < u, w < |Vu|. (3) 

The above setting arises naturally for inertial particles 
in an incompressible turbulent flow u(£, x) where parti- 
cle's drag by the flow is describable by the linear Stokes 
friction (see [9J for details). The Newton law reads 



q = v, tv + v = u(t, q[t] 



(4) 



At r — > the particle follows the flow, Vq = u(t, q[t] 
For leading order correction in r put v = Vq in tv, 



v(t) wm - r[dtu + (u ■ V)tt] 



(5) 



where vq = dtu + (u ■ V)u is the substantial derivative 
along the trajectory of the fluid particle and the RHS is 
evaluated at x = q(t). In this order the particle velocity 
is determined by its location in space. One may introduce 
the velocity field v(t,x) = u — r[dfU + (u ■ V)it] that 
determines particles trajectories in space as in Eq. ©. 
The small correction v — u gives v a finite compressibility 



= -rV • [(it • V)«] = tV 2 P ^ 0. 



(6) 



Above we took divergence of the Navier-Stokes equations 
dtu + (it • V)tt = — Vp + vV 2 u + f where p is the pres- 
sure, v is the viscosity and the force / is assumed to be 
either incompressible or concentrated on the boundaries. 
One can use only the potential part of the correction, 
v = u + rVp, as it becomes clear from the formulas 
below. Formation of a singular density is characteristic 
for particles and it holds already for a centrifuge: par- 
ticles shift a bit after each rotation eventually reaching 
the boundary to form a singular density. Incompressible 
turbulence can be thought of as a collection of vortices- 
centrifuges where particles tend to accumulate on the 
boundaries between the vortices, cf. 

Unless v is degenerate usrb is non-smooth. The de- 
generacy was quantified in Q, based on a Green-Kubo 
type formula for the space-averaged sums of backward- 
A~ and forward-in-time A^~ Lyapunov exponents [3, |8| , 
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dt / w(0,x)w[t,q(t,x)]dx. (7) 



At weak compressibility one may substitute q(t, x) above 
by the trajectories of incompressible flow X(t, x) denned 
by d t X(t, x) = u[t, X(t, x)] and X(0, x) = x. We define 
the "microcanonical" correlation function by the usual 

(w(0)w(t)) = w(Q,x)w[t,X(t,x)]dx. (8) 



By incompressibility (w(0)w(t)) = (w(0)w(—t)) so that 

/oo 
(w(0)w(t))dt/2 - -E(0)/2, 
-oo 

where E(0) is the spectrum of w[t, X(t,r)] at zero fre- 
quency, E(0) > 0. Assuming the sums of Lyapunov expo- 
nents are equal for all r except for a set of zero volume 
(e.g. this holds for systems obeying the SRB theorem 
043)) and remembering the sums determine the rates of 
growth of density 0] , we have for the solutions to Eq. ((J} , 

hm lln^^=-hm 1 In ^ 
Moot n(0,r) t->oo t n[0,q{0\t,r)] 2 

where q(0\t,r) is the trajectory passing via r at time t. 
For smooth n(t — 0), the above implies that, in a non- 
degenerate case of E(0) > 0, at large times, n(t) tends 
to zero at all r except for a set of zero volume, while be- 
coming infinite on q(t, r). This rules out smooth usrb 
and shows a strange attractor is formed for Eqs. ©-© 
generally. It is immediate to infer the Kaplan- Yorke codi- 
mension A of this attractor. At weak compressibility we 



have Aj < but 



d-l 



A; 



-A? > 0, where we as- 



sume that the Lyapunov exponents A° of u do not vanish 
identically (this implies A^ < by £ A° = 0). Thus 



A = ]>>/A d «£:(0)/(2|A o d |). 



(9) 



Applying the above to particles in turbulence, we find 

/oo 
(<K0Mt))dt = -r 2 £ (O)/2 < 0, 
-oo 

where £^,(0) > is the spectrum of </> = V 2 p and we 
used that there is no degeneracy. Thus inertial particles 
in turbulence asymptote a fractal set with the Kaplan- 
Yorke codimension A = r 2 ^(0)/(2|A§|) where Aj_ < is 
the third Lyapunov exponent of turbulence [Tol . Il2| . 

Phase space density serves to find averages of dynam- 
ical variables - functions on the phase space, so only in- 
tegrals of usrb with smooth test functions are observ- 
able. Thus "weak solutions" to the continuity equation - 
possibly singular nsRB{r) satisfying J drnsRB(r)v(r) ■ 
^/( r ) — for all smooth f(r) (coinciding with V • 
[n-SRBv] = if usrb is smooth) describe physically 
meaningful states. Our main assumption is that such 
solutions are obtainable as a long-time limit T — > oo of 
evolution of a unit initial density set at t = —T: 

n S R B (r)= lim e p{ ~ T ^ r \ p(t,r)= f w[q(t', r)]dt'. (10) 

T^oo J Q 

The factor p(t, r) describes changes in an infinitesimal 
volume V(t) located at t = near r, as p(t,r) = 
ln[V(t)/V(0)] where t can be both positive and negative. 
In limT^oo j dre p{ - T < r \r)v{r)-Vf(r) first one takes the 



integral over r and then the limit T — >• oo. Then Eq. (|T0|) 
defines a weak solution usrb if for any smooth /, 

lim J/(T)=0, If(T)=f dre p{ - T > r) v(r)Vf(r). (11) 

T-j-oo J 
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Changing variables r = q(T,x), using Vi[q(T,x)] = 
Vj(x)\7jqi(T,x), see @, and integrating by parts, one 
finds If(T) = -(w(0)/(T)). We will assume that corre- 
lations decay (w(0)/(T)) = (w(0))(/(T)) = at T -> oo 
and cj[q(t, r)] has a finite correlation time r c . Then 
Eq. (|10p defines a weak solution to the continuity equa- 
tion. This obeys the SRB theorem [2j on the equality of 
time and phase-space averages for dynamical variables /, 

lim (1/T) / f[q(t,x)]dt= f f(r)n SRB (r)dr, (12) 

provided the time-average is independent of x for al- 
most every x. The proof is based on the identity 
limr-s-oo J f[q(T,x)]dx = J f{r)n SRB (r)dr. The den- 
sity defined by Eq. (jlOp is an exponent of a sum of a 
large, asymptotically infinite, number ~ T/t c of uncor- 
rclatcd random variables. This implies hsrb is singular 
pointwise. However, its correlation ip(r) at r ^ is finite, 

v(r) = (n SRB (0)nsRB(r)}= lim / e /»(-r,o)+ P (-r,r)\ (13) 

T— >oo\ I 

where the angular brackets designate spatial averaging. 
The cumulant expansion theorem [l|, [l3[ , applied to the 
volume conservation {exp[p(t, r)]) = 1, gives the "sum 
rule" Y^=x {p n {^i r )) c / n ' = 0' where the subscript c 
stands for the cumulant (giving to lowest order in com- 
pressibility — (p) = (y0 2 ) c /2). Applying the expansion to 
Eq. (|13|) and using the sum rule, one finds only the mixed 
terms, containing both p(t, r) and p(t, 0) do not cancel, 

oo 

ln*>(r)= lim £[<[pi W) c -<tf> c -(^> c ]M (14) 

n=2 

where pi = p(—T, 0) and pi = p(—T, r). The above leads 
to a series representation for the attractor's correlation 
dimension, cf. 14]. In the limit of small compressibility, 

r° 

^(r) = e 9 « g(r)= dtidt^wpCfa, 0)]w[X(i 2) r)]).(15) 

t/ — oo 

Above we put X(t, r) instead of q(t, r) and omitted the 
subscript c since for incompressible flow (w["X(t, r)]) = 
(w(r)) and (w(r)) — J V • txir = 0. The condition of 
applicability of Eqs. (|T5j) is the negligibility of n > 2 
terms in Eq. (| 14[) . Proceeding for (71(1*1)71^2). .n(r/v)) = 
lim r _ i . 00 (exp[X)p(-T,r i )]) as for ip(r), 

\n( n SRB(ri)n S RB{r2)--nsRB{rN)) = ^.g(r 4 - r 3 ).(16) 

i>j 

Thus TigRB has log-normal statistics completely deter- 
mined by a single structure function g(r). The latter 
behaves universally at small r. It diverges at r = be- 
cause X(t, 0) and X(t, r) do not separate and spend in- 
finite time together. Thus (n 2 SRB ) — oo reflecting tisrb 
is singular and compressibility is a singular perturba- 
tion. At small but finite r the trajectories separate expo- 
nentially with the Lyapunov exponent |A°| of the time- 
reversed flow u. The trajectories stay together during 



a time that diverges logarithmically at r — > 0. To sin- 
gle out this divergence fix a scale r, where below 
r] both v and uj are smooth, and consider r such that 
t* = ln[r»/r]/| A° 3> t c . Separating the interval of inte- 
gration over ti in Eq. (fT5|) into (— oo, — i*) and (— i*,0), 

/■— t* /*oo 

ff(r) = / dti/ dt 2 (a;[X(t 1 ,0)]a;[X(t 2 ,r)])+^ J B(0).(17) 

J —oo J —oo 

The first term on the RHS is estimated as g(r») while the 
divergence at small r is due to the last term: 

cp(r) « (r*/r) 2A . (18) 

Thus the correlation codimension is twice the Kaplan- 
Yorke one. Note 77 is defined up to a factor of order unity 
but this is of no effect for the answer by A <C 1. The 
statistics is isotropic at small r independently of isotropy 
of v. The origin of this is that the separation vector 
between two infinitesimally close trajectories grows with 
an exponent independent of its initial direction. 

Applying the above to the inertial particles, the statis- 
tics is completely determined by the structure function 

S(r) = f dMi 2 (#i,X(t 1 ,0)Mi 2 ,X(t 2 ,r)]>, 

J — OO 

involving only the properties of turbulence. For density, 

{n S R B (0)n S R B (r)) = exp[r 2 6»]. (19) 

To analyze the domain of validity of Eq. (jT9")l con- 
sider n = 3 term in Eq. (fT^)l . It is of the type 
r 3 J dt x dt 2 dU{<l>[h, X(ti, 0)]# 2 , X(t 2 , 0)M* 2 , X(* 2 , r)]) c . 
Its evaluation within Kolmogorov's theory (KT) |15| 
shows it is smaller than r 2 S(r) by the factor of the 
Stokes number St = t| A§ | , the product of r and a 
typical value of the velocity gradient. This leads to 
the validity condition St <C 1 of Eq. (fH)]) . Due to the 
deviations from KT, higher n terms in Eq. (|14[) contain 
additional powers of the Reynolds number Re and the 
true condition is StRe" <C 1 with a > 0. One expects 
a <C 1 so in practical situations St <C 1 ensures Eq. (fH?j) . 
The same condition is expected to satisfy the additional 
demand that the underlying Eq. ([5]) holds. For particles, 
X] ~ (^/lA^I) 1 / 2 is the Kolmogorov scale of turbulence 
where v the fluid viscosity [la ] . Neglecting intermittency, 
T 2 S(rj) ~ St 2 so that g{rj) <C 1 and one can neglect the 
first term in Eq. (Tl7j) for t* = ln(?,/r)/|A§|, 

(nsRB(0)nsRB(r)) = [rj/r) , r < ??. (20) 

The above was obtained in [lol ] . There are no significant 
density fluctuations, {tisrb (fi)nsRB (**)) ~ («Si?.s) 2 = 1 
already at r <C 7/, where A In (r//r) <C 1. Turbulence ef- 
fect on the particles depends on the observer's resolution: 
at A In (77/7- ) > 1 segregation holds, while at larger r - 
mixing. This is how the assumption of mixing can work 
effectively for particles on a multifractal. Segregation 
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can also bring physical effects regardless the observer's 
resolution by influencing the rate of collisions [10|]. 

At St 1 density inhomogeneities are absent in the 
inertial ran ge r ) < r < L, where L is the pumping scale of 



turbulence [15| . Then Eq. (|20|) is a complete description. 
In contrast, at St ~ 1 the inertial range inhomogeneities 
are important [l8[ and our Eq. p^|) . extended to hold 
asymptotically at St ~ 1, gives a unique access to the 
inhomogeneities. In KT S(r) depends only on r and the 
mean energy dissipation e so \n(p KT (r) fts C(t/t^ t ) 2 , 
where C is dimensionless constant and t^ t = e _1 / 3 r 2 ' 3 
is the Kolmogorov time-scale of turbulent fluctuations at 
scale r. The KT prediction applies rigourously to the 
Kraichnan model of v decorrelated in time [l2l . Hint, 
However for turbulence the numerical simulations [18| 



show lny(r) 



oc r 



-10/3 _ This is at moderate Re where 



KT is expected to work. Noting S(r) ~ T 2 cy([p(x + r) — 
p(x)] 2 }, where r r is the relevant time-scale, we suggest 
the difference has the same origin as the deviations of 
the pressure scaling from KT |19| . 

The absence of fluctuations beyond Aln^/r) > 1 is 
the central property of weakly dissipative dynamical sys- 
tems. Consider the density ni coarse-grained over the 
volume Vi of a ball with radius I in d dimensions, 



(r) = mi(r)/Vi, mi 



n SR B(r')dr', (21) 



|r'-r|<( 



The mass m; is equal to the mass contained in the dy- 
namical image of the ball at time —t which is an ellip- 
soid with largest axis growing as I exp[|A°|i] (fluctuations 
are negligible by weak compressibility). For — t = = 
\a{R/l)/\\° d \ where R < ?y, but Aln(7i/i?) < 1, this mass 
is spread over a region where the density self-averages and 
it equals just the volume of that region, Vi exp[p(— t*, r)], 



ni(r) = exp[ / o(-<*,r)]. 



(22) 



regime: density fluctuates at scale I due to mass con- 
densing from a volume with size smaller than r\ over 
which the mass is effectively distributed uniformly. Us- 
ing -(p 2 (-t)) n /2 = (p(-t))=tT.XT and (v/R) 2A w 1, 



-(p 2 M))c/2 = ( P (-t)) =tJ2K ^d (v/R) 2 ^ 
one finds the lognormal statistics (rtf) = (rj/l)^^ x \ ci 
Eq. HU). To find the fractal dimensions D(a) [HIH]], 



D(a)i 



:limln(m" 1 n S RB)/[{a 



l)lnZ] (23) 



consider (n" tisrb) 
f:£u[t,q{t,r)}di\). 
p(-t*,r) w p{-U H 
dent 



averaging \exp 



= lim T _ >00 (exp[ap(-i st , r) + 
For ^> r c , using 

r c ,r) one can make indepen- 



This is a fundamental formula of the weakly dissipative 



ap(-U,r) + j_ T * u(q(t, r) 

(exp [ap(— t, r)])(exp J_^u(q(t,r) ) = (nf), where we 

used volume conservation. This gives D(a) — d — Aa: 
fractal dimensions are close to d and the attractor is 
almost space-filling (for a 3> 1 and high moments the 
lognormal approximation is not valid generally). The 
general result of [2lj for D{a) in the Kraichnan model 
in d = 2 reduces to our result at small compressibility. 

Our results give testable predictions for particles in 
real turbulence with no need for assumptions about the 
turbulence structure. The fractal structure, with dimen- 
sions defined by a single parameter A, exists below r cor 
and is formed in a characteristic time ln(?7/r COT .) /| A3 1 ~ 
A3| 2 /St 2 i?0(O). An everyday application is the distri- 
bution of water droplets in warm clouds, important for 
such problems as radiation or rain formation in clouds 
[Toj . The analysis can be extended to include gravity 
and the case of light particles and to describe correla- 
tions between densities of different size (r) particles. 
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